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ABSTRACT 

The character of star formation is intimately related to the supersonic magne¬ 
tohydrodynamic (MHD) turbulent dynamics of the molecular clouds in which stars 
form. A significant amount of the turbulent energy dissipates in low velocity shocks. 
Fast and slow MHD shocks differ in how they compress and heat the molecular gas, 
and so their radiative signatures reveal distinct physical conditions. 

We use a two-fluid model to compare one-dimensional fast and slow MHD shocks 
propagating at low speeds (a few kms~^). Fast shocks are magnetically driven, forcing 
ion species to stream through the neutral gas ahead of the shock front. This magnetic 
precursor heats the gas sufficiently to create a large, warm transition zone where all 
the fluid variables smoothly change in the shock front. In contrast, slow shocks are 
driven by gas pressure, and neutral species collide with ion species in a thin hot slab 
that closely resembles an ordinary gas dynamic shock. 

We consider shocks at velocities Vg = 2-4kms~^ and preshock Hydrogen nuclei 
densities nn = 10^-10^ cm“^. We include a simple oxygen chemistry and cooling by 
CO, H 2 and H 2 O. CO rotational lines above J = 6 ^ 5 are more strongly excited in 
slow shocks. These slow shock signatures may have already been observed in infrared 
dark clouds in the Milky Way. 
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1 INTRODUCTION 

Understanding the internal environment of the giant molec¬ 
ular clouds (GMCs) in which stars form is a neces¬ 
sary precursor to addressing the s tar formation rate and 
stellar initial mass fun ction (e.g., iBergin fc Tafallal l2007l : 
iMcKee fc Qstril^l2007ll . Larg e non-therma l linewidths ob¬ 
serve d in molecular lines fe.g.. [bar^ nl ll98ll : ISolomon et al.l 
[1983) have often been attributed to turbulent motions. Fur¬ 
thermore, the kinetic energy associated with these motions 
is generally found to be on the order of the gravitational po¬ 
tential energy, indicating the importance of turbulence as a 
dynamical component of molecular clouds. In addition, mag¬ 
netic fields in GMCs give Alfven velo cities on the order of 
the observed velocity dispersions (e. g. JCrutcher et'^ll993l : 
ICrutcherlll999l : rCrutcher et al.ll2010ll . This suggests that the 
turbulence in GMCs is magnetohydrodynamic (MHD) in na¬ 
ture. 

The physics underlying the MHD turbulence of molec¬ 
ular clouds is intimately connected with properties of 


* E-mail: andrew.lehmann@mq.edu.au 


star formation J^a c Low fc Klesseiil 200j) . For exam¬ 
ple. iFederrath fc KlessenI ( 2012l l and iFederrath fc KlessenI 
(120131) use three dimensional MHD simulations to analyse 
different modes of turbulence by comparing compressive 
driving to solenoidal driving. They showed that the star for¬ 
mation rate and efficiency are both sensitive to these driving 
modes and Mach number variations. Hence it is desirable to 
discover observable distinctions between different modes of 
turbulence. 

Supe rsonic turbulence dissip ates via shock waves and 
vortices llPetv fc Fafgarona 2000li. In s i mulat ions of com¬ 
pressible MHD turbulence IStone et al.l (Il998l i found that 
shock waves dissipated 50% of the turbulent energy in 
a strong magnetic field model and 32% in a weak field 
model. The heating of compressed gas in a thin post¬ 
shock region uniquely drives chemistry and radiative cool¬ 
ing. Thus the radiative characteristics of turbulent dissipa¬ 
tion will be strongly shap ed by cooling in shocks. Further¬ 
more, [SmhLeLky (l2000bh found that weak shocks in a large 
range of velocities were responsible for the majority of dissi¬ 
pation in simulations of decaying MHD turbulence. In con¬ 
trast, a small range of stronger shocks dissipated turbulence 
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that was being driven dSmith et al.ll2000^ i. Hence radiative 
signatures of shocks could be used to distinguish between 
these two scenarios. 

The observational signatures of the low-velocity shocks 
that dominate dissipati on of MHD turbu lence have only re¬ 
cently been considered. IPon et ^ ] (l20l3i considered C-type 
fast MHD shocks travelling at speeds of 2-3kms“^ per¬ 
pendicular to the magnetic field and computed the abun¬ 
dances and emission of H 2 , CO and H 2 O. By comparing 
CO rotational line emission from these shocks to those 
produced in photodissociation regions (PDRs), they found 
that fast s hocks dominate the em ission in transitions above 
J = 5 4. iLesaffre et al.l (l2013tl take a statistical approach 

by computing observational diagnostics due to a distribu¬ 
tion of C- and J-type fast, perpendicular, MHD shocks at 
velocities ranging from 3-40 km s~^. They use these shocks 
to explain the radiation from a turbulent wake formed by a 
galaxy collision in Stephan’s Quintet. 

Anomalously bright CO lines above J = 5 —>■ 4 
have been recen tly observed towards Milky Way molecu - 
lar clouds (e.g., |Pon et al.| |20l4 |2015|; iLarson et ahl |2015l~) 
and from warm molecular gas in external g alaxies (e.g., 
iKamenetzkv et ^l2012l : iPellegrini et al.ll2013l i. These stud¬ 
ies all conclude that PDR models are unable to reproduce 
the bright high-J CO lines and all suggest MHD shock waves 
as the heating mechanis m. The shock models used or re¬ 
ferred to in these studies jFlower fc Pineau Des Foretsll20f^ : 
|Pon et ^ l2012f) and other studies of shocks in interstel- 
lar cl ouds (e.g. jDraine_&Kut41l98^; iHollenbach fc McKed 
1 19891 : IChapuian fc Wardld 2009 ') all consider fast MHD 
shocks. MHD fluids can, however, support three kinds of 
shocks: fast, intermediate and slow. Unfortunately no study 
has identified which kinds of shocks would be produced by 
MHD turbulence. One of the goals of our work is to moti¬ 
vate the classification of MHD shocks in turbulent molecular 
clouds. 

Ideal MHD assumes a fully ionized gas in which the 
magnetic field is frozen. Molecular clouds are in fact only 
weakly ionized and the magnetic field acts on the neutral 
fluid via the ionized fluid. A two-fluid description is there- 
for e more appropriate to molecu lar cl oud studies. Ear l y work 
by iLithwick fc Goldreichl ll200ll ') and iLazarian et al.l (l2004l l 
highlighted the importance of two-fluid effects on the turbu¬ 
lent cascade. Recently it has become feasible to run high res¬ 
olution t hree dimensional tw o-fluid MHD simulations. For 
example, iMever et al.l ll2014l ~l use simulations to show that 
observations of linewidth differences between line emission 
from neutral and ion spec ies can be accounted f or by two- 
fluid effects. Furthermore, iBurkhart et al.l (l2015l i use simu¬ 
lations to show that the Alfvenic modes do not necessarily 
dissipate at the ambip olar diffusion sc ale. They also confirm 
the an alytic results of iBalsai^ (Il996tl and iTillev fc Balsaral 
(l201lll that in molecular cloud conditions, the fast and 
Alfven modes are strongly damped by ion-neutral collisions 
leaving the slow modes to propagate with little damping. 
Thus nonlinear steepening into slow shocks might be ex¬ 
pected to preferentially occur in molecular clouds. In addi¬ 
tion, slow shocks reach far higher peak temperatures than 
fast shocks for two reasons. Firstly, the heating timescale in 
fast shocks is determined by the long ion-neutral collision 
timescale while in slow shocks is determined by the short 
neutral-ion collision timescale. As the cooling timescale of 


the gas lies between these two, the gas being overrun by fast 
shocks remains at low temperatures whereas it quickly heats 
up within slow shocks. Secondly, as we show in Section[2]the 
gas in fast shocks necessarily loses some of its kinetic energy 
to strengthening the magnetic field while in slow shocks it 
does not. Thus slow shocks have more energy available to 
heat the gas. 

In this paper, we solve the steady, plane parallel two- 
fluid MHD equations to model shocks that propagate at any 
angle to the magnetic field. In Section [2] we elucidate some 
of the basic differences of fast and slow MHD shocks. We de¬ 
scribe our computational scheme in Section [31 In Section [3] 
we compare fast and slow shocks in the low-velocity regime 
with molecular cloud conditions. We discuss their radia¬ 
tive characteristics and implications for interpreting emis¬ 
sion from turbulent molecular gas. Finally, in Section [5] we 
show how these signatures might be used to interpret obser¬ 
vations. 


2 THEORY 

Cosmic rays streaming through the interstellar medium 
weakly ionize molecular clouds, generating an ion fluid that 
interpenetrates the neutral particles. A two-fluid MHD de¬ 
scription is developed in Section [2T] with the ion fluid cou¬ 
pled to the magnetic field via the Lorentz force and to the 
neutral fluid via the collisional force, so that 


J X B 

c 


apipn (Vi — v„) 


where a is the rate coefficient for elastic ion-neutral scatter¬ 
ing and p and v are the density and velocity with subscripts 
i and n referring to ion and neutral fluids respectively. By 
considering the state of these fluids far away from the shock 
front, we illustrate how the various families of MHD shocks 
come about and highlight the different effects they have on 
the ambient magnetic field. 

A simplihed chemical model, in which two body re¬ 
actions and photodissociation affect the abundances of 
coolants, is presented in Section 12.21 Finally, the details of 
how those coolants radiate are outlined in Section [2]3] Only 
rotational excitations of molecular species are considered be¬ 
cause vibrational excitations do not become signihcant until 
temperatures reach > 1000 K, which are not achieved in the 
low velocity shocks computed here. 


2.1 Two-Fluid MHD 

Conservation equations 

While t urbulence is an inhere ntly three-dimensional prob¬ 
lem fe.g. [Burkhart et al.ll20l4 and references therin), shock 
waves are extremely thin structure s within a turbulent sys¬ 
tem. Hence we follow IPrainel (Il986l l by considering a station¬ 
ary plane-parallel shock travelling in the 2 direction with the 
magnetic field initially lying in the x-z plane. The governing 
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equations of the ion fluid are 


^ ipiVi.) - 0 , 

d f B^B^\ _ 

^ I ^ 1 — Olpipn [Vix Vnx) 

d (Bl\ _ 

d^i'^ ~ 


( 1 ) 

( 2 ) 

(3) 


The neutral fluid equations are 


— (pnVnz) — 0, (4) 

(pnVnzVnx) ~ Otpipn (Uix Vnx') , (5) 

^ H” -Pn) — Otpipn {Viz Vnz) ; (6) 


where Pn is the neutral pressure, F is the heating func¬ 
tion, A is the cooling function, and the internal energy 
u = Pn/ (7 — 1) for adiabatic index 7 . Finally, the electro¬ 
magnetic equations give 


j i^VizBx VixBz) — 0 , 

dz 

( 8 ) 


(9) 


Solving for the internal structure of intermediate shocks re¬ 
quires equations analogous to Equations m to for the 
y-direction. We consider only fast and slow shocks here by 
ignoring this case. Far away from the shock front there is no 
velocity difference between the fluids. This means there is 
no frictional heating in those regions, and so the shocks sat¬ 
isfy the isothermal one-fluid jump conditions. The types of 
shocks allowed are therefore determined by the MHD signal 
speeds. 


Shock Families 


When considering small perturbations of density, pressure, 
the velocity field, and the magnetic field around time in¬ 
dependent averages the equations of MHD allow for three 
linear wave modes: the fast, intermediate and slow waves. 
These waves travel with phase velocities 


/ = + ^\j{v\ + cif - 4:v\ci cos2 9^ 

^ \/(«i + cos 2 


i = va cos 9, 


2 , 2 

'Ca + Ca 


1/2 


1/2 


where va = Bly/inp is the Alfven velocity, Cs = ^/ksT/pm 
is the isothermal sound speed with Boltzmann constant fcs 
and mean mass per particle prn ~ {7/S)mH, and 9 is the 
angle between the magnetic field and the direction of prop¬ 
agation under consideration. These speeds are plotted as 
functions of 9 in Figure [T] for the case va > c^ (the relevant 
case in molecular clouds). 

In the frame of reference comoving with a shock front, 
the preshock fluid travels toward the shock front at a 
shock velocity Vs greater than one of the wave speeds. The 
three wave speeds demarcate four regions of fluid velocities 



Figure 1. The phase velocities of linear MHD wave modes versus 
the angle between the magnetic field and direction of propagation 
of those modes for va > Cs- The wave speeds delineate the regions 
marked 1 to 4. 


marked 1 to 4 in Figure [T] Inside of the shock front the fluid 
must transition down across a wave speed, which allows for 
six kinds of shock waves collected into three families: fast, 
intermediate and slow MHD shocks. Fast shocks are those 
that cross the fast wave speed only (1 —>■ 2 ), intermediate 
shocks cross the intermediate wave speed (1 —^ 3, 1 —^ 4, 
2—^3 and 2 —>■ 4), and slow shocks cross the slow wave 
speed only (3—^4). 

Far ahead and far behind a shock where there is no 
friction between the fluids—and hence v„ = Vi—one can 
show from Equations © to ([9]) that the product 



is conserved across any shock. This can be used to emphasise 
some basic differences between the different kinds of shocks. 
Suppose, for simplicity that > 0 in the preshock medium. 
Recall that a fast shock crosses the fast wave speed only, so 
that the term in parentheses remains positive but is reduced 
across the shock. Thus B^ must increase to compensate. In 
intermediate shocks, the velocity crosses the intermediate 
wave speed and so the term in parentheses switches sign 
across the shock. This implies that Bx must switch sign 
also. Finally, in slow shocks, the velocity crosses the slow 
wave speed only, so that the term in parentheses is negative 
and becomes further negative in the postshock medium. This 
means that Bx must decrease to compensate. These three ef¬ 
fects on the magnetic field direction are shown schematically 
in Figure [2] 

The switch in sign of Bx in the intermediate shock is 
due to a rotation of the magnetic field within the shock 
front. As the field rotates out of the x-z plane, intermediate 
shocks require equations analogous to Equations m to 
for the 3 /-direction. This case is ignored here because it is 
unclear whether stead y state intermediate shocks are phys i- 
cally admissable (e.g., IWuIIiQStI : iFalle fc Komissaro^ 1200 l| j . 
Furthermore, the 1 —>■ 3 intermediate shock will resemble a 
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Fast Shock Intermediate Shock Slow Shock 



Figure 2. The effect on magnetic field orientation of the three 
classes of MHD shock waves. Fast shocks (left) increase the angle 
between the field and shock normal, intermediate shocks (middle) 
reverse the sign of the angle and slow shocks (right) decrease it. 


fast shock as it crosses the fast speed, the heating in a 2 4 

shock front will be dominated by a hydrodynamic jump as 
it crosses the sound speed—see Section [3]— and thus resem¬ 
ble a slow shock, the 2—^3 shock will only weakly heat the 
gas as it resembles a rotational discontinuity, and finally the 
1—^4 shock will resemble a fast shock followed by a weak 
2^-3 shock followed by a slow sh ock, so our models s hould 
also roughly capture its structure dKennel et al.l[l989l) . 

As the spacing of field lines is proportional to the mag¬ 
netic field strength, one can see in Figure [2] that the field 
strength increases across fast shocks and decreases across 
slow shocks. This means that some of the kinetic energy of 
a fast shock is converted into magnetic field energy. Hence, 
for slow shocks at the same velocity one expects there to be 
more energy available to heat the gas. 


Governing Differential Equations 

With some manipulation. Equations m to reduce to 
three ordinary differential equations determining changes in 
the X component of the magnetic field and the neutral fluid 
temperature: 

dBa: vlapioBl ( Vi^ - v„ 


vlB. 


VnzViz 


dz 

dTn 

dz Pnouf (uL - ) ' 


Tu (7 - 1) 


( 10 ) 

( 11 ) 


— Vs ) (r — A) — apipn {Viz — Vnz) VnzVt 


with the velocities 

v\Bz 


VsBl 


{Bx — Bxo) : 


'^ix — 


VizBx VsBxO 

~^z ’ 


Vnz = y (^/3 ± 'A- 4r^ , 

VsBxqBx -f Bz i^BxVnx BzVnz^ 


Viz = 


the preshock Alfven velocity 


VA = 




Bo 




and 


^ _l_ kaTno 1 Bio — Bl 


ks 


p.mVi 


PmVl ' 2 vl Bl 

■Tn 


where a subscript 0 denote preshock values. 

Equations and (HU are not complete without spec¬ 
ifying how energy leaves the gas via A and enters via F. The 
collisions between the ion and neutral fluids generate fric¬ 
tional heating at a rate 

Ff = apnPi (Vi - v„)^ . 

per unit volume. In addition, we include cosmic-ray heating, 
which is only important at temperatures of 10-30 K. 

The shock heated gas cools by radiating away heat en¬ 
ergy stored in the rotational and vibrational modes of its 
constituents. These modes are collisionally excited, so that 
the cooling function must be a function of the densities of the 
coolants, the density of the colliding particles and the tem¬ 
perature of the gas. The abundances of coolants can change 
due to chemical reactions, and so we present a simple model 
for the oxygen chemistry occuring in the shock heated gas 
in the next section before specifying how that gas radiates 
in Section 


2.2 Chemistry 


Radiative cooling depends on the the abundances of the 
coolants and hence chemical reactions influence the cool¬ 
ing of the gas. The domina nt coolants in molecular c louds 
are CO, H 2 and H 2 O (e.g., iNeufeld fc Kaufmanlll993l) and 
the abundance of the latter is determined in this study 
by t he set of eight r e action s listed in Table [T] adopted 
from IWagner fc Grafj (Il987l ). Previous models of the gas 
phase chemistry in interstellar environments—using more 
than 100 reactions— were found to be dom inated by a small 
set of reactions ( e.K., llgIesias fc Silklll978l ). The reactions of 
IWagner fc Grafll are the subset of these dominant reactions 
that control the abundance of H 2 O. 

The change in the particle density of a neutral species 
M through the shock is given by 

y {n{M)vnz) = Sm 

where Sm is the rate at which M is created or destroyed. 
X) Sm over all neutral species must be zero in order for 
Equation to remain true. The low densities of molecular 
clouds mean that only two body chemical reactions and pho¬ 
todissociation processes need to be considered. Two-body 
reaction rates between species A with particle density n{A) 
and species B with particle density n{B) take the form 

kAB{T)n{A)n{B) cm~® s~^ 


where kAB (T) is a rate coefficient . The UMIST Datab ase for 
Astrochemistry 2012 (RATE12) llMcElrov et al.ir2013l ) gives 
this coefficient in the form 


feAfl(r) = a 


( T 

= rv ( - 


V300 


exp 


(-S 


3 -1 
cm s 


where a, /3 and 7 are constants. The values of these pa¬ 
rameters for the reactions used are shown in Table [T] These 
rates are verified for the temperature range 200 — 2500 K, 
so there is some uncertainty when extrapolating down to 10 
K. In fact, for the reaction 


T = 


O -k OH O 2 + H 

the rate coefhcient parameters in the RATE12 database give 
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Table 1. Reaction rate coefficient parameters. 


No. 

Reaction 

a 

/3 

7 

1 

0 - 1 - H 2 OH + H 

3.14 X 10-13 

2.70 

3150.0 

2 

OH -I- H 0 -1- H 2 

6.99 X 10-11 

2.80 

1950.0 

3 

OH-I-H 2 -)-H20-I-H 

2.05 X 10-13 

1.52 

1736.0 

4 

H 2 O-1-H ^ OH-1-H 2 

1.59 X 10-11 

1.20 

9610.0 

5 

OH-I-OHH 2 O-1 -0 

1.65 X 10-13 

1.14 

50.0 

6 

H 2 O-1 -0 ^ OH-1-OH 

1.85 X 10-11 

0.95 

8571.0 

7 

0 -1- OH O 2 -f H 

4.33 X 10-11 

-0.5 

30.0 

8 

O 2 -1-H ^ 0 -1-OH 

2.61 X 10-1° 

0 

8156.0 

Notes. 

The parameters for 

reaction 7 are 

taken 

from 

IWaerner Sz Graffl JlQSTI^ instead of RATE12. 



a rate that diverges in the 

10-100 K range. 

For this rate 


the parameters are taken from IWagner fc Grafll (Il987l l. In 
Fignreini the rates are plotted against temperature. It can be 
seen that most of the rates only “tnrn on” at temperatnres 
T > 60 K, which is why the weak fast shocks considered 
here only negligibly affect the molecular abundances. 

Photodissociation in molecular clouds is caused by ul¬ 
traviolet radiation generated by secondary electrons in cos¬ 
mic ray ionization events (?). The rates of photodissociation 
events per volume of species A with particle density n{A) 
take the form 


R { ... —3 —1 

PA- - 'n{A) cm s 

1 — cj 

where C,cr is the cosmic ray ionization rate, oj is the albedo 
of the dust grains found in molecular clouds and pa is an effi¬ 
ciency constant. The ionization rate t^cR is set to 10”^^ s~^, 
the albedo w is set to 0.6, and the values of pa fo r the 
reactions considered—taken from iGredel et ^ (Il989ll — are 
shown in Tabled 

With these two rates and using the static, planar as¬ 
sumptions the abundance of species M relative to the hy¬ 
drogen nuclei density xm = n(M)/nH changes as 

where no is the preshock H-nuclei density and the sums are 
taken over reactions that either produce or destroy M. 

The end result of this chemistry is to add five more dif¬ 
ferential equations—coupled to each other through densities 
and to Equations and m through temperature—to 
follow the abund ances of H, O, OH, O 2 and H 2 O. We follow 
|Pon et ^ 1 i|2oI 3) in using an initial H 2 O abundance of 10 
O abundance of 5.45 x 10“"^ and a C abundance of 1.4 x 10“"^ 
which is assumed to be entirely locked up in CO. The CO 
abundance is assumed to be constant throughout the shock, 
because it has a dissociation temperature higher than any 
temperature reached in these weak shocks. The initial H, 
OH and O 2 abundances are set to 10“"^, and 

respectively. H 2 abundance is then computed through the 
shock using xh 2 = (1/2) (1 — xh)- 


2.3 Cooling 

The cooling function of iNeufeld fc KaufmanI lll993ll is 
adopted for this model. This uses an escape probability to 
account for the effects of reabsorption by the surrounding 



Figure 3. Reaction rate coefficients as a function of temperature. 


Table 2. Photodissociation rate efficiency. 


Reaction 

PA 

H 2 O -1- Photon OH -I- H 

971 

OH -1- Photon —> 0 -1- H 

509 

O 2 + Photon ^0 + 0 

751 


Notes . Photodissociation rate efficiencies taken from lOredel et ^ 
lll989ll . 


media on the rotational level populations. Collisions with 
H 2 are the only excitations considered as the abundance of 
H 2 in molecular clouds is orders of magnitude above the 
next most abundant molecular species. The power radiated 
by CO, H 2 O and H 2 in a wide range of conditions is ex¬ 
pressed in terms of a rate coeffecient Lm defined such that 
the power radiated per unit volume by species M is 

Am = n(M)n(H2)I/M. 


Lm is a function of H 2 density, temperature and an optical 
depth parameter N{M). N{M) is a correction factor that 
accounts for reabsorption of radiation by the same molecule 
in the surrounding gas. It depends on the geometry of the 
system in question, the density of M and the local velocity 
gradient. The expression for a plane-parallel slab of thickness 
d and characteristic velocity difference An, 


fV(M) 


n(M)d 

9An 


( 12 ) 


is used here. Lm is then expressed with a four parameter 
analytic fit to its density dependence 


_L = _L + ^(H2) / n(H 2 ) y L _ ni/ 2 Lo \ 

Lm Lo Llte Lo \ n-^/2 J V Llte J 


where Lq is the low density limit of the cooling rate coef¬ 
ficient, Llte is the luminosity per molecule with the level 
population in local thermodynamic equilibrium and ni /2 is 
the H 2 density at which the cooling rate coefficient is half of 
Lq. Lo is a function of temperature only while Llte, ni/2 
and a are functions of temperature and N{M). 
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Preliminary isothermal shocks were computed in order 
to gain an estimate on the values of d and An, and these 
values checked for self consistency against the results when 
cooling was included. The average value of n{M) through 
the shock was used. H 2 is optically thin in all interesting 
astrophysical conditions, so Lh 2 does not depend on Equa¬ 
tion (II 2 II . 

The values of Lp, Llte, n,i /2 and a for CO, H 2 O and 
H 2 are taken from lNeufeld et ^ l|l99 5ll in the temp e rature 
range T = 10-100 K, and from Neufeld fc Kaufmad (Il993l l 
in the range T = 100-2000 K. We combine the parameters 
in the T = 10-100 K range using an ortho-to-para H 2 O ratio 
of 3:1, which was already assumed in the T = 100-2000 K 
range. The values of some of these parameters at T = 100 
K are inconsistent between the two sources, so we take the 
mean value. 


3 NUMERICAL INTEGRATION 


The core of this problem are the first order ordinary differ¬ 
ential equations which take the form: 


dBx 

dz 

dz 


fi{z,B,,T) 

f 2 {z,B,,T) 


for fi and /2 defined by Equations cni) and m- They must 
be solved simultaneously with the abundance derivatives 
which are coupled to these two through their dependence on 
temperature. When supplied with initial conditions, these 
equations can be integrated to give Bx, T and the abun¬ 
dances as functions of 2 . The open source PYTHON mod¬ 
ule sakits.ode^ was used, which solves initial value prob¬ 
lems for ODEs using variable-order, variable-step, multistep 
methods. 

The initial conditions are stationary states—in that 
their 2 derivatives are zero—so integrating from these points 
changes nothing. A perturbation must be added to these ini¬ 
tial conditions in the form 


Bx — BxO + SBx€ ^ 
T = To + STe^'^^ 


before integrating. The effect of the perturbation can be 
understood by linearising the differential equations da and 
m- It is illuminating to look at the isothermal case—where 
7 = 1—so that only the B field derivative remains. In this 
case we get the eigenvalue 

. ^ apio jvs - f) {vl - s^) 

Vs v\ {vl - ci) 

where /, i and s are the fast, intermediate and slow sig¬ 
nal speeds defined in Section [ 2 T] Xx determines whether 
the perturbation of Bx grows or decays. By replacing the 
preshock variables with their postshock counterparts, the 
eigenvalue can also be used to explore the region of solution 
space near postshock states. 

For fast shocks, the preshock state is in region 1 of Fig¬ 
ure [ 1 ] so that Vs > f > i > Cs > s. This means Xx is 


^ https://github.com/bmcage/odes 


positive so that Bx grows. This describes an unstable sta¬ 
tionary point, where any perturbation away from the initial 
condition grows. The postshock state is in region 2 so that 
f>Vf>i>Cs>s where Vf is the final velocity. This 
means Xx is negative so that this state is a stable stationary 
point. In Section [2.II we noted that the fast shock increases 
Bx, hence a positive perturbation of the initial Bx is all that 
is required to finish at the fast postshock state. 

For slow shocks, the preshock state is in region 3 of Fig¬ 
ure [ 1 ] so that f > i > Vs > Cs > s in the supersonic case. 
This means Xx is negative and this state is a stable station¬ 
ary point. Hence there is no way to leave the supersonic 
slow preshock state in a continuous fashion like in C-type 
fast shocks. The slow solution is further complicated by a 
singularity in the equations when crossing the sound speed 
{vnz —>■ Cs), as a manipulation of Equation d3 gives 

dVnz _ (7 1) (r A) apfipi {Viz Vfiz) Vjiz 

dz Pn{cl-vlz) 


We cross this sonic point by inserting a gas dynamic jump 
in the neutrals determined by the hydrodynamic jump con¬ 
ditions: 


V2 


Vl 

Tl 

Ti 


7-1 

7 -I- 1 
1 -b 


+ 


7-b 1 M2 


M^ (7 - 1 ) -b 2 
M2 (7 -b 1) 


(13) 

(14) 


where M^ = p\v\MP\ and V 2 and T 2 are the neutral 2 
velocity and temperature immediately after the gas dynamic 
jump. The slow postshock state lies in region 4 of Figure [T] 
so that f>i>Cs>s>Vf and Xx is negative. This means 
this state is a stable stationary point, and so jumping across 
the sound speed (via Equation (1131) 1 will allow the solution 
to smoothly settle onto the slow postshock state. 

For slow shocks with a subsonic preshock state, f > i > 
Cs > Vs > s so that Ax is positive and the stationary state is 
an unstable. This means a perturbation that reduces Bx will 
grow smoothly until the solution reaches the slow postshock 
state. Such a shock is a C-type slow MHD shock, and will 
be ignored here because it requires high sound speeds and 
therefore high temperatures that are not obviously relevant 
for molecular cloud studies. 


4 RESULTS 

Here we compare a set of shock models with parameters 
shown in Table [3] The turbulent cascade of energy en¬ 
sures the dissipation i s dom inated by low-velocity shocks, so 
we follow |Pon et akl ll2012|j in looking at velocities around 
3kms“'^. In Figure 0] the Alfven velocity is plotted against 
the magnetic field orientation 6 with dashed lines of con¬ 
stant intermediate speed i = va cos 0 overlaid. The magnetic 
field strength on the second vertical axis is computed for a 
preshock density no = 10® cm~®. The intermediate speed 
is the upper bound of possible slow shock velocities—region 
3 in Figure [T]— and so the dashed lines in Figure 0] trace 
the minimum field strength required to get slow shocks with 
peak temperatures shown. We thus chose 0 low enough to re¬ 
tain reasonable field strength values. At 9 — 30°, a 4kms“^ 
slow shock—which will heat the gas to ~ 857 K immedi¬ 
ately after the neutral jump—requires an Alfven velocity of 
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Figure 4. Alfven velocity versus magnetic field orientation, with 
dashed lines of constant intermediate speed va cos 9 that al¬ 
low slow MHD shocks to reach the peak temperatures shown. 
The magnetic field strength on the second vertical axis assumes 
preshock density no = 10® cm~®. 


~ 4.6kms“^, which is fixed for all the slow shocks at differ¬ 
ent preshock densities. This gives magnetic field strengths 
in the range Bq = 25-253 pG. Finally, N was initially cho¬ 
sen for preliminary isothermal shocks, and then recomputed 
after each shock model for consistency. Table [3] shows the 
self consistent values of N(CO) and N(H20). 


Structural characteristics 

In Figures[S]and[n]we compare fast and slow shock profiles of 
ion and neutral velocity, density, temperature, cooling rate 
and abundances for Va = 3kms“^. The fast shock propa¬ 
gates at 89.9° to a 10 pG magnetic field, while the slow 
shock propagates at 30° to an 80 pG field. The preshock 
density no = 10® cm“® for both. 

In the velocity profiles—upper panels of Figures [S] and 
[E] —the neutral and ion velocities in the shock propagation 
direction are the dashed and dotted lines respectively. These 
plots are in the frame of reference of the shock wave, and 
so velocities below the shock velocity represent fluid flowing 
ahead of the shock front in the lab frame, in which the shock 
wave travels to the left. Therefore, in the fast shock the ions 
stream ahead, imparting some of their momentum to the 
neutrals until both fluids are moving at the same speed in 
the postshock medium (on the right side of each plot). This 
process is set by the long ion-neutral collision timescale and 
therefore happens smoothly over a large distance, and there 
is no viscous jump as seen in hydrodynamic shocks. In con¬ 
trast, for the slow shock it is the neutrals that flow ahead 
of the ions that then get accelerated to the neutral veloc¬ 
ity over a small distance—set by the fast neutral-ion colli¬ 
sion timescale—giving slow shocks a much thinner structure 
than the fast shocks. The cooling timescale of the gas lies 
between the two collision timescales, and so slow shocks will 
reach higher peak temperatures than fast shocks at the same 
velocity. 


Table 3. Shock Parameters. 


Parameters 

Vs 

km/s 

logN(CO) 
log ^cm“^ 

logN(H20) 

(km/s)-®) 

Fast Shocks, 6 = 89.9° 

no = 

10® cm ® 

2.0 

15.9 

11.7 

Bo = 

3 uG 

2.5 

15.7 

11.8 



3.0 

15.6 

11.9 



3.5 

15.5 

11.9 



4.0 

15.5 

11.9 

no = 

10® cm“® 

2.0 

15.8 

12.7 

Bo = 

10 uG 

2.5 

15.8 

12.7 



3.0 

15.8 

12.6 



3.5 

15.8 

12.6 



4.0 

15.5 

12.3 

no = 

lO"' cm-® 

2.0 

15.9 

12.8 

Bo = 

32 g-G 

2.5 

15.8 

12.7 



3.0 

15.8 

12.8 



3.5 

15.7 

12.6 



4.0 

15.7 

12.6 

Slow Shocks, 6 = 30° 

no = 

10® cm ® 

2.0 

14.5 

11.0 

Bo = 

25 uG 

2.5 

14.5 

11.0 



3.0 

14.5 

11.0 



3.5 

14.5 

11.4 



4.0 

14.5 

12.3 

no = 

10® cm”® 

2.0 

14.5 

11.0 

Bo = 

80 g-G 

2.5 

14.5 

11.0 



3.0 

14.5 

11.5 



3.5 

14.5 

12.6 



4.0 

14.5 

13.2 

no = 

lO"' cm-® 

2.0 

15.6 

12.6 

Bo = 

253 g.G 

2.5 

15.7 

12.7 



3.0 

15.1 

13.6 



3.5 

15.0 

14.6 



4.0 

15.0 

15.1 


The density profiles—solid black lines in the upper 
panels—similarly show large differences between the shocks. 
The fast shock weakly compresses the gas by a factor of ~ 7 
in a simple, smooth manner. The slow shock has a complex 
density profile that ends with a compression ratio of ~ 300 
and must be understood in conjunction with the tempera¬ 
ture and cooling (see middle panel) taking place. It starts 
with the neutral jump, governed by Equations GSl) and d, 
compressing the gas by a factor of ~ 4 and heating it to 
~ 500 K. Combining Equations ([3]) and we get 


d 

dz 



= 0 


where we identify each term as effective pressure terms: ram 
pressure Pram — PriV^z, gas pressure Pgaa = Pn = nksT 
and magnetic pressure Pmag = /Sn. As the neutrals lead 
the ions in the slow shock. Equation implies that the 
magnetic pressure must drop. The velocity drop causes the 
ram pressure to drop, and so gas pressure must increase, 
which is shown in the third panel of Figure [fi] Efficient line 
emission from CO and H 2 causes the temperature to quickly 
decrease, so to compensate the density must increase. This 
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Figure 5. Profiles for a = 3kms“^ fast MHD shock propa- 
gating through a medium with initial number density no = 10® 
cm“® at 89.9° to a 10 |xG magnetic field. The first panel shows 
the velocity and density profiles, the second panel shows the tem¬ 
perature and cooling rate profiles, the third panel shows the gas, 
ram and magnetic pressure profiles, and the fourth panel shows 
the abundances of selected oxygen molecules. 


occurs smoothly in the ~ 0.1 x 10^® cm after the initial 
jump, until the density reaches a plateau at a compression 
ratio of ~ 300. The density finally settles as the neutral and 
ion velocities equalise. 

While the velocity and density structural differences are 
large, the observational implications are sensitive to the tem¬ 
perature and cooling rate profiles shown in the second panels 


Figure 6. Same as Figure [5l but for a Us = 3kms“® slow shock 
propagating through a medium with initial number density no = 
10® cm~® at 30° to an 80 pG magnetic field. 

of Figures [H and [6l These panels show that the peak temper¬ 
atures differ by an order of magnitude with the slow shock 
peaking at ~ 487 K and the fast shock peaking at ~ 50 K. 
At this temperature, the fast shock causes CO emission only 
slightly above the background level. The slow shock is hot 
enough to emit significantly in fines of CO and H 2 , and also 
in H 2 O above the CO background emission. The fourth pan¬ 
els of Figures [ 5 ] and [6] also show that the slow shock drives 
more chemistry than the fast shock due to its high tempera¬ 
ture. The strong H 2 O emission in the slow shock is partially 
due to this increase in its abundance. 


Temperature (K) 
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Following the pressure terms—in the third panels of 
Figures [5] and [6] —through the shocks provides an intuitive 
understanding of the fundamental differences between fast 
and slow shocks. In the fast shock, the magnetic pressure 
(dotted lines) is higher in the postshock region than in the 
preshock region. As the gas pressure (solid lines) is neglible 
everywhere in this shock it is understood to be magneti¬ 
cally driven, where the high Pmag region pushes into the 
low Pmag region. The ion fluid is strongly coupled to the 
magnetic field, and so the ions are pushed forward ahead of 
the neutrals. This is why the ions lead the neutrals in the 
velocity profile. This situation is reversed in the slow shock, 
where a high gas pressure region pushes the neutrals through 
the ions. In this case, the ion coupling to the magnetic field 
deforms it and increases the separation between field lines, 
which reduces the held strength. 

Comparison of fluxes 

The cooling discussed in Section 12.31 can be used to search 
for observational differences between the kinds of shocks. 
Integrating the cooling rate—middle panels of Figures [5] and 
[5] — through the shocks gives the hux emitted in the direction 
of the shock normal. In Figure [T] we compare the CO, H 2 
and H 2 O fluxes of all the fast and slow shocks in Tabled 

The CO flux is similar between both kinds and across 
the velocity range. The H 2 hux is much more temperature 
sensitive and is therefore stronger in the slow shocks by 5 
orders of magnitude at the lowest velocity and 1 order of 
magnitude at the highest. The H 2 O hux is also stronger in 
the slow shock at all velocities. 

Figure [7] shows that for slow shocks the CO and H 2 
huxes are within an order of magnitude of each other at all 
velocities. For fast shocks, however, the CO hux is always 
stronger by more than an order of magnitude. This feature 
holds at initial densities no = 10^ and 10"^ cm“®, though the 
magnitudes of the huxes are lower and higher respectively 
by factors of 10 than the no = 10^ cm“^ case. This suggests 
that combining observations of hux from both molecules is 
a strong indicator of the kind of shock being observed. 


Rotational lines 

Figure [7] shows similar levels of CO hux for both fast and 
slow shocks at all the velocities considered here. However, 
the strength of the individual rotational lines that make 
up this hux is not expected to be the same given the 
higher temperatures that slow shocks reach. Here we esti¬ 
mate those those line strengths using the non local thermo- 
dynamic equilibrium (LT E) radiative transfer code RADEX 
dvan der Tak et al.ll2007l l. 

For a given radiating molecule, RADEX requires as in¬ 
put the density of H 2 as the collisional partner, the column 
density of the radiating molecule, the temperature and the 
linewidth. We consider 4 molecules with H 2 as the colli¬ 
sional partner: ^^CO, ^®CO, ortho-H20 and para-H20. We 
use a ^^CO to '^^CO ra tio of 61:1, w e ll wit hin the range of 
measurements given bv iMilam et al.l (|2005h . and an ortho- 
to-para H 2 O ratio of 3:1 as assumed in Section [2.31 for the 
cooling function. 

We use RADEX in slab mode, and use density weighted 



2.0 2.5 3.0 3.5 4.0 

Vg (km/s) 


Figure 7. Energy flux contributions by different coolants in 
shocks of initial density no = 10^ cm“^ (top)? '^0 = 10^ cm~^ 
(middle) and no = 10^ cm~^ (bottom), and shock velocities 
Vs = 2-4kms~^. Blue circles refer to fast shocks and red tri¬ 
angles refer to slow shocks. Solid, dotted and dashed line types 
indicate whether the coolant is CO, H 2 or H 2 O respectively. 
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Figure 8. Estimates of integrated line fluxes from rotational lines of (red), (green), ortho-H20 (blue) and para-H20 (black) 

for the fast (top) and slow (bottom) shocks shown in Figures [5] and |6] 


averages of the required inputs over appropriate slab defi¬ 
nitions of the computed shocks. For slow shocks, we define 
multiple slabs in order to account for the complex tempera¬ 
ture and density structure. While RADEX accounts for op¬ 
tical depth effects within each slab and outputs the optical 
depth for each rotational transition, the emission from a slab 
may be reabsorbed by the other slabs it has to pass to reach 
the observer. Hence estimates of line strengths from lines of 
high optical depth—such as the CO lines below J = 5 —>■ 4— 
may not be reliable. Differences in the optical depth of par¬ 
ticular lines in fast and sl ow shocks may have a n observa¬ 
tional effect. For example. iBurkhart et al.l (l2013l ) show that 
the slope of the spatial power spectrum derived from syn¬ 
thetic observations of J = 2 1 in simulations of 

MHD turbulence is sensitive to the optical depth of this 
line. Furthermore, this slope could distinguish between sub- 
and super-Alfvenic turbulence. However, to relate the opti¬ 
cal depth of a line emitting from dense postshock regions to 
that over a line of sight through a cloud requires knowledge 
of the spatial distribution of shock waves in the cloud, which 
is beyond the scope of this paper. Finally, we normalise the 
line strengths so that emission from each molecule from all 
its lines is equal to the fluxes given in Figured 

In Figure [S] we plot the line strengths as computed by 
RADEX for the fast and slow shocks shown in Figures[5]and 
[6l As expected, the higher temperature of the slow shock 
results in strong excitation of the high-J CO transitions of 
both isotopes. There are also excited H 2 O lines in the slow 
shock that are negligible in the fast shock. 

Any line of sight through a cloud will intersect multiple 
shocks propagating at different velocities, so we look for fea¬ 
tures in the spectra that hold across the velocity range. In 
FigureOwe plot the line integrated fluxes for selected high- 
J lines of CO across the velocity and density ranges for fast 
and slow shocks. In all slow shocks, the dynamic range of 
these high-J lines of CO is much lower than in fast shocks. 
For example, the line ratio CO J = 5 ^ 4/10 —>■ 9 for a 
fast shock is always greater than 30 times the ratio from the 
slow shock of the same velocity and density. 

As line ratios don’t change if the line strengths are re¬ 
duced by a constant factor, these distinguishing features will 
remain even though the shocks propagate in various direc¬ 
tions with respect to the line of sight. Furthermore, any line 


of sight through a turbulent cloud will cross multiple shock 
fronts. Emission from optically thin lines will be a simple 
addition of emission from each shock, retaining line ratio 
characteristics which therefore could be strong indicators of 
shock type. 

The high temperatures reached within slow shocks allow 
H 2 cooling to become comparable to cooling by CO, as can 
be seen in Figure[7l The H 2 molecule has no dipole moment, 
and so this cooling is due to weak quadrupole emission. We 
can estimate the low lying pure rotational lines (v = 0—^0) 
of this emission using Figure 1 of iBurton et al.l (1 19921 ). which 
assumes the column of radiating H 2 is in LTE. For pure ro¬ 
tational lines S(0)-S(3), the line integrated fluxes emitted 
normal to the plane of the slow shock in Figure [6] are in 
the range ~ erg/cm^/s. These lines a strongly 

suppressed in gas below 100 K, and so the fast shock in 
Figure [5] produces negligible S(0)-S(3) line emission. This 
estimate shows that the pure rotational lines of H 2 could 
be an important diagnostic of shock type. However, the as¬ 
sumption of LTE doesn’t always hold and so to produce a 
more accurate prediction for H 2 rovibrational line emission 
the level populations would have to be computed in parallel 
with the shock flow variables and reaction rat e equations, as 
is don e for instance in the fast shock model of iGusdorf et al.l 

(l2008l) . 


5 DISCUSSION 


We have integrated the two-fluid MHD equations to ob¬ 
tain one dimensional time-independent fast and slow shock 
wave solutions. While two-fluid fast MHD shocks have been 
well studied (e.g. iDrain 3 l 19861 : iFlower fc Pineau Des Forit3 
Il998l) , two-fluid slow shocks have not been considered in 
molecular cloud conditions. Thus we use a simplified model 
in order to highlight the qualitative differences between 
the two kinds of shocks. Simplifications include restrict¬ 
ing the chemical network to a small subset of reactions 
that influence the abundance of H 2 O, as this is an im¬ 
portant coolant in molecular gas. Gas-phase oxygen chem- 
istry is well containe d in the restricted network used here 
lllglesias fc_S^ 197^ and so a richer network—such as that 
used in ( Glover et al. llioTa) —would not strongly affect H 2 O 
production. We have not modeled how a variety of differ- 
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Figure 9. Integrated line intensities of CO rotational lines for 
fast (blue circles) and slow (red triangles) shocks with preshock 
densities no = 10^ cm“® (top)i no = 10® cm“® (middle) and 
no = 10“^ cm“® (bottom), and shock velocities Vs = 2 —4kms~®. 
The solid lines connect rotational lines of the same transition from 
J = 5 -h 4 to J = 10 -h 9. 


ent initial abundance s—which varies th roughout a turbu¬ 
lent molecular cloud llGlover et alJlioi^ —could affect the 
shock chemistry or cooling profiles. However, this will not 
change the general differences in structure between fast and 
slow shocks which is determined by whether the shock is 
driven by magnetic or gas pressure, respectively. Rotational 
line emission from H 2 could provide more observational di- 
agnostics and has b een computed in fast shock models by 
iLesaffre et al.l (|201.ll ) . In order to predict this emission accu¬ 
rately the level populations of the rovibrational states of H 2 
have to be computed in parallel with the shock flow variables 
and reaction rate equations. 

This work shows that fast and slow MHD shocks are 
structurally and observationally distinct. In Section [2.11 we 
noted that fast shocks increase the angle between the mag¬ 
netic held and direction of propagation whereas slow shocks 
decrease this angle. In principle this effect could be observed 
in studies of the polarized thermal emission revealing the ge- 
ometry of magnetic fields in molec ular clouds. The study by 
[Planck Collaboration et al.i (l2015l ) has shown that in high 
column density filaments the magnetic field tend s to be 
perpendicular to the filament. IPadoan et al.l (1200 ill explain 
these filaments as dense postshock regions resulting from the 
collision of supersonic turbulent flows. This result precludes 
the possibility that the filam ents contain the fast shocks 
most commonly mode led (e.g. IPlower fc Pineau Des Forets! 
I 2 OI 0 I : |Pon et ^ 12012 !) in which the magnetic field is paral¬ 
lel to the shock front. Higher spatial resolution work of this 
kind could either detect the magnetic held bending across 
shock fronts or rule out the shock formation scenario of fil¬ 
aments. In Section |4] we showed that observations of high-J 
CO lines—above J = 6 —>■ 5—could distinguish between 
fast and slow MHD shocks. Multiple transitions should si¬ 
multaneously be observed to disentangle the possible shock 
models through the use of line ratios. In addition, veloc¬ 
ity information from linewidths could be used to constrain 
the shock velocity of the model. Hence, useful observations 
would require spectral resolutions better than 1 km/s at the 
high-J CO lines (CO J = 6 —^5toJ = 10—>-9 lines lie in 
the frequency range 691.47-1151.99 GHz). In addition, the 
spatial resolution must be sufficient to avoid gas warmed by 
protostellar outflows or stellar winds. By way of example, 
we consider a recent observation that satisfies these criteria 
and show that slow shock signatures may be present. 

|Pon et ^ (l2015l) (hereafter P15) present observations 
of CO J = 8 —>-7, 9— ^8 and 10 —>■ 9 taken with the 
Herschel Space Observatory, towards four starless clumps 
within Galactic infrared dark clouds (IRDCs). These clumps 
were chosen because they lacked massive embedded proto¬ 
stars, avoiding confusion from outflows which could also cre¬ 
ate a warm gas component. The authors detect CO J = 8 —>■ 
7 and 9 —>■ 8 towards three of their clumps—named Cl, FI 
and F2—and give an upper limit for the 10 —>■ 9 line. They 
compare these observations to PDR models at densities of 
10'* and 10®cm“®, typical of IRDCs, and interstellar radi¬ 
ation fields of 1 and ~ 3 Habing. All of the PDR models 
underpredict the CO J = 9 ^ 8 line, so the authors sug¬ 
gest that the dissipation of turbulence in low velocity shocks 
could account for these lines. Here we consider how fast and 
slow shocks similar to those modeled in this paper could 
account for these observations. 

To convert the integrated line intensity (f Tdv) of P15 
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Table 4. CO high-J line luminosities for IRDCs 


ijPon et al ][20ll) and selected shock models 


Source 

CO J = 8 -)• 7 

^10-14 erg/s/cm^^ 

CO J = 9 ->■ 8/8 ->■ 7 

CO J = 10 9/8 7 

Cl 

1.26 ±0.06 

0.65 ±0.07 

< 2.10 ±0.24 

FI 

1.48 ± 0.09 

0.94 ±0.10 

< 2.20 ±0.29 

F2 

1.19 ±0.09 

0.77 ±0.10 

< 2.69 ±0.36 

Slow Shock Models 

(A) 

1.66 

0.51 

0.27 

(B) 

1.75 

0.80 

0.58 

Fast Shock Models 

(A) 

1.19 

0.16 

0.03 

(B) 

4.10 

0.25 

0.06 


Notes. The values used for the J = 10 ^ 9 line from IRDCs are the suggested upper limits. Slow shocks A and B have shock velocities 
Vs = 2kms~^ and 3.5kms“^ respectively. Fast shocks A and B have shock velocities Vs = 3.5kms“^ and 4.0kms“^ respectively. All 

four shock models have preshock density np = 10“^ cm~^. 


to flux we use the formula 

P=^jTdv 

where k is the Boltzmann constant, is the beam area, u 
is the frequency of the transition nnder consideration and 
c is the speed of light. The CO J = 8—>-7, 9—>-8 and 
10 9 transitions have rest frequencies of 921.80, 1036.91 

and 1151.99 GHz and half power beam widths (HPBWs) of 
23, 20 and 19 arcseconds respectively. We use 

^ 

to compute the beam area. Finally, we use the average of 
the detected line intensities for comparison (Column 8 of 
Table 3 in Ponl5). For this set of values the line fluxes are 
shown in Table S) The second column shows the CO J = 
8 —>■ 7 line flux, while the third and fourth columns show the 
CO J = 9 8/8 —>■ 7 and CO J = 10 —>■ 9/8 —>■ 7 line ratios 

respectively. The CO J = 10 9 line was not detected in 

any of the clumps so we use the upper limits adopted in 
Ponl5. This means that the CO J = 10 —^ 9/8 —>■ 7 line 
ratio is an upper limit. 

We also list in Table|4]the predicted values from selected 
slow and fast shocks by multiplying the RADEX intensities 
(Figure |9]) by f2/47r for the appropriate beam areas. The 
predicted value of the CO J = 8 —^ 7 integrated intensity 
assumes the shock front faces the observer and fills the beam. 
We chose the models that give the closest CO J = 8 ^ 7 
line fluxes (Slow and Fast shock A in Table |4]) as well as the 
shocks that give the closest CO J = 9—>-8/8—>-7 line ratios 
(Slow and Fast shock B in Table |4]). 

Slow shocks A and B have shock velocities Vs = 2 km s“^ 
and 3.5 kms”*^ respectively. Fast shocks A and B have shock 
velocities Vs = 3.5kms~ and 4.0 kms”*^ respectively. All 
four shock models have preshock density no = 10'* cm“® 
which agrees well with typical densities of IRDCs. Fast shock 
A is the only shock that doesn’t overpredict the observed in¬ 
tegrated intensities, but the predicted integrated intensities 
of the other three models could be reduced to match the 
observations if the shock doesn’t fill the beam. Hence the 
observed line flux of CO J — 8 ^ 7 from these IRDCs can 
be explained by either fast or slow shocks. 

The slow shocks generally fit the CO J = 9 8/8 7 

line ratios better than the fast shocks, with the line ratio 


from slow shock B (0.80) very close to the average of the 
ratio for all three clumps (~ 0.78). The predicted line ra¬ 
tio from fast shock B (0.25) is the largest predicted from 
all the fast shocks modeled, which still underpredicts the 
observed values by ~ 3 — 4 times. Hence the high tempera¬ 
tures produced in slow shocks are necessary to explain both 
the CO J = 8 ^ 7 and 9 —>■ 8 emission from these clumps. 
Combined with the inability of PDR models to explain these 
observations, we therefore suggest that a slow shock inter¬ 
pretation is favoured by the models in this paper. 


6 CONCLUSIONS 

The one dimensional time independent two-fluid MHD equa¬ 
tions were numerically integrated to compare the structure 
of low-velocity fast and slow shocks in molecular clouds. Our 
simplified model includes the effects of the major coolants 
found in molecular clouds and follows the abundances of 
chemicals affecting the production of H 2 O in a simple chem¬ 
ical network. The solutions highlight important differences 
between fast and slow MHD shocks in molecular clouds. 
These shocks show strong differences in their velocity and 
density structure because of the different driving pressures 
behind the shock fronts. Fast shocks are driven by magnetic 
pressure while slow shocks are driven by gas pressure. This 
means that the thickness of fast shocks is set by the long ion- 
neutral collision timescale, wheareas the thickness of slow 
shocks is set by the short neutral-ion collision timescale. 
The cooling timescale of the gas lies between these two and 
so peak temperatures in slow shocks are far higher than fast 
shocks of the same shock speed. 

'We showed that fast and slow shocks are observation- 
ally distinct and provided some example diagnostics. For in¬ 
stance, low lying pure rotational lines of H 2 contribute neg¬ 
ligibly to the cooling in fast shocks, whereas they produce 
significant radiation from the warm gas in slow shocks. The 
non-LTE radiative transfer code RADEX was used to esti¬ 
mate line strengths of rotational transitions of *^CO, *®CO, 
ortho-H 2 0 and para-H 2 0. The higher temperatures of slow 
shocks excite the high-J transitions of ^^CO more than in 
fast shocks. Line ratios near these transitions show strong 
differences between fast and slow shocks across the velocity 
range and therefore may be strong indicators of shock type. 










Fast and slow MHD shock waves 13 


Anomalously strong high-J CO lines have been observed 
in nearby infrared dark clouds (|Pon et al.ll2015ll . The line 
ratios from these observations closely match slow shock pre¬ 
dictions and are poorly fit by any of the fast shocks modeled 
here. 

This suggests that simulations of MHD turbulence 
could gain observational predictions if the statistics of shock 
types were recorded. If the mixture of shock families is found 
to be sensitive to turbulence parameters—such as the driv¬ 
ing mode, the kind of feedback included, Mach number vari¬ 
ations and self-gravity—then shock signatures become ob¬ 
servational probes of the turbulence. Combined with tracers 
of star formation like young stellar objects, shock signatures 
could then shed light on the influence of supersonic MHD 
turbulence on the character of star formation. 
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